3.473 \(\int \frac{(a^2+2 a b x^2+b^2 x^4)^3}{x^{21}} \, dx\)

Optimal. Leaf size=84 \[ \frac{b^3 \left (a+b x^2\right )^7}{1680 a^4 x^{14}}-\frac{b^2 \left (a+b x^2\right )^7}{240 a^3 x^{16}}+\frac{b \left (a+b x^2\right )^7}{60 a^2 x^{18}}-\frac{\left (a+b x^2\right )^7}{20 a x^{20}} \]

[Out]

-(a + b*x^2)^7/(20*a*x^20) + (b*(a + b*x^2)^7)/(60*a^2*x^18) - (b^2*(a + b*x^2)^7)/(240*a^3*x^16) + (b^3*(a +
b*x^2)^7)/(1680*a^4*x^14)

________________________________________________________________________________________

Rubi [A]  time = 0.056049, antiderivative size = 84, normalized size of antiderivative = 1., number of steps used = 6, number of rules used = 4, integrand size = 24, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.167, Rules used = {28, 266, 45, 37} \[ \frac{b^3 \left (a+b x^2\right )^7}{1680 a^4 x^{14}}-\frac{b^2 \left (a+b x^2\right )^7}{240 a^3 x^{16}}+\frac{b \left (a+b x^2\right )^7}{60 a^2 x^{18}}-\frac{\left (a+b x^2\right )^7}{20 a x^{20}} \]

Antiderivative was successfully verified.

[In]

Int[(a^2 + 2*a*b*x^2 + b^2*x^4)^3/x^21,x]

[Out]

-(a + b*x^2)^7/(20*a*x^20) + (b*(a + b*x^2)^7)/(60*a^2*x^18) - (b^2*(a + b*x^2)^7)/(240*a^3*x^16) + (b^3*(a +
b*x^2)^7)/(1680*a^4*x^14)

Rule 28

Int[(u_.)*((a_) + (c_.)*(x_)^(n2_.) + (b_.)*(x_)^(n_))^(p_.), x_Symbol] :> Dist[1/c^p, Int[u*(b/2 + c*x^n)^(2*
p), x], x] /; FreeQ[{a, b, c, n}, x] && EqQ[n2, 2*n] && EqQ[b^2 - 4*a*c, 0] && IntegerQ[p]

Rule 266

Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Dist[1/n, Subst[Int[x^(Simplify[(m + 1)/n] - 1)*(a
+ b*x)^p, x], x, x^n], x] /; FreeQ[{a, b, m, n, p}, x] && IntegerQ[Simplify[(m + 1)/n]]

Rule 45

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[((a + b*x)^(m + 1)*(c + d*x)^(n + 1
))/((b*c - a*d)*(m + 1)), x] - Dist[(d*Simplify[m + n + 2])/((b*c - a*d)*(m + 1)), Int[(a + b*x)^Simplify[m +
1]*(c + d*x)^n, x], x] /; FreeQ[{a, b, c, d, m, n}, x] && NeQ[b*c - a*d, 0] && ILtQ[Simplify[m + n + 2], 0] &&
NeQ[m, -1] &&  !(LtQ[m, -1] && LtQ[n, -1] && (EqQ[a, 0] || (NeQ[c, 0] && LtQ[m - n, 0] && IntegerQ[n]))) && (
SumSimplerQ[m, 1] ||  !SumSimplerQ[n, 1])

Rule 37

Int[((a_.) + (b_.)*(x_))^(m_.)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[((a + b*x)^(m + 1)*(c + d*x)^(n +
1))/((b*c - a*d)*(m + 1)), x] /; FreeQ[{a, b, c, d, m, n}, x] && NeQ[b*c - a*d, 0] && EqQ[m + n + 2, 0] && NeQ
[m, -1]

Rubi steps

\begin{align*} \int \frac{\left (a^2+2 a b x^2+b^2 x^4\right )^3}{x^{21}} \, dx &=\frac{\int \frac{\left (a b+b^2 x^2\right )^6}{x^{21}} \, dx}{b^6}\\ &=\frac{\operatorname{Subst}\left (\int \frac{\left (a b+b^2 x\right )^6}{x^{11}} \, dx,x,x^2\right )}{2 b^6}\\ &=-\frac{\left (a+b x^2\right )^7}{20 a x^{20}}-\frac{3 \operatorname{Subst}\left (\int \frac{\left (a b+b^2 x\right )^6}{x^{10}} \, dx,x,x^2\right )}{20 a b^5}\\ &=-\frac{\left (a+b x^2\right )^7}{20 a x^{20}}+\frac{b \left (a+b x^2\right )^7}{60 a^2 x^{18}}+\frac{\operatorname{Subst}\left (\int \frac{\left (a b+b^2 x\right )^6}{x^9} \, dx,x,x^2\right )}{30 a^2 b^4}\\ &=-\frac{\left (a+b x^2\right )^7}{20 a x^{20}}+\frac{b \left (a+b x^2\right )^7}{60 a^2 x^{18}}-\frac{b^2 \left (a+b x^2\right )^7}{240 a^3 x^{16}}-\frac{\operatorname{Subst}\left (\int \frac{\left (a b+b^2 x\right )^6}{x^8} \, dx,x,x^2\right )}{240 a^3 b^3}\\ &=-\frac{\left (a+b x^2\right )^7}{20 a x^{20}}+\frac{b \left (a+b x^2\right )^7}{60 a^2 x^{18}}-\frac{b^2 \left (a+b x^2\right )^7}{240 a^3 x^{16}}+\frac{b^3 \left (a+b x^2\right )^7}{1680 a^4 x^{14}}\\ \end{align*}

Mathematica [A]  time = 0.0078792, size = 82, normalized size = 0.98 \[ -\frac{15 a^4 b^2}{16 x^{16}}-\frac{10 a^3 b^3}{7 x^{14}}-\frac{5 a^2 b^4}{4 x^{12}}-\frac{a^5 b}{3 x^{18}}-\frac{a^6}{20 x^{20}}-\frac{3 a b^5}{5 x^{10}}-\frac{b^6}{8 x^8} \]

Antiderivative was successfully verified.

[In]

Integrate[(a^2 + 2*a*b*x^2 + b^2*x^4)^3/x^21,x]

[Out]

-a^6/(20*x^20) - (a^5*b)/(3*x^18) - (15*a^4*b^2)/(16*x^16) - (10*a^3*b^3)/(7*x^14) - (5*a^2*b^4)/(4*x^12) - (3
*a*b^5)/(5*x^10) - b^6/(8*x^8)

________________________________________________________________________________________

Maple [A]  time = 0.048, size = 69, normalized size = 0.8 \begin{align*} -{\frac{{a}^{5}b}{3\,{x}^{18}}}-{\frac{15\,{a}^{4}{b}^{2}}{16\,{x}^{16}}}-{\frac{5\,{a}^{2}{b}^{4}}{4\,{x}^{12}}}-{\frac{3\,a{b}^{5}}{5\,{x}^{10}}}-{\frac{{a}^{6}}{20\,{x}^{20}}}-{\frac{10\,{a}^{3}{b}^{3}}{7\,{x}^{14}}}-{\frac{{b}^{6}}{8\,{x}^{8}}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((b^2*x^4+2*a*b*x^2+a^2)^3/x^21,x)

[Out]

-1/3*a^5*b/x^18-15/16*a^4*b^2/x^16-5/4*a^2*b^4/x^12-3/5*a*b^5/x^10-1/20*a^6/x^20-10/7*a^3*b^3/x^14-1/8*b^6/x^8

________________________________________________________________________________________

Maxima [A]  time = 1.0085, size = 95, normalized size = 1.13 \begin{align*} -\frac{210 \, b^{6} x^{12} + 1008 \, a b^{5} x^{10} + 2100 \, a^{2} b^{4} x^{8} + 2400 \, a^{3} b^{3} x^{6} + 1575 \, a^{4} b^{2} x^{4} + 560 \, a^{5} b x^{2} + 84 \, a^{6}}{1680 \, x^{20}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((b^2*x^4+2*a*b*x^2+a^2)^3/x^21,x, algorithm="maxima")

[Out]

-1/1680*(210*b^6*x^12 + 1008*a*b^5*x^10 + 2100*a^2*b^4*x^8 + 2400*a^3*b^3*x^6 + 1575*a^4*b^2*x^4 + 560*a^5*b*x
^2 + 84*a^6)/x^20

________________________________________________________________________________________

Fricas [A]  time = 1.69188, size = 174, normalized size = 2.07 \begin{align*} -\frac{210 \, b^{6} x^{12} + 1008 \, a b^{5} x^{10} + 2100 \, a^{2} b^{4} x^{8} + 2400 \, a^{3} b^{3} x^{6} + 1575 \, a^{4} b^{2} x^{4} + 560 \, a^{5} b x^{2} + 84 \, a^{6}}{1680 \, x^{20}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((b^2*x^4+2*a*b*x^2+a^2)^3/x^21,x, algorithm="fricas")

[Out]

-1/1680*(210*b^6*x^12 + 1008*a*b^5*x^10 + 2100*a^2*b^4*x^8 + 2400*a^3*b^3*x^6 + 1575*a^4*b^2*x^4 + 560*a^5*b*x
^2 + 84*a^6)/x^20

________________________________________________________________________________________

Sympy [A]  time = 0.902203, size = 75, normalized size = 0.89 \begin{align*} - \frac{84 a^{6} + 560 a^{5} b x^{2} + 1575 a^{4} b^{2} x^{4} + 2400 a^{3} b^{3} x^{6} + 2100 a^{2} b^{4} x^{8} + 1008 a b^{5} x^{10} + 210 b^{6} x^{12}}{1680 x^{20}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((b**2*x**4+2*a*b*x**2+a**2)**3/x**21,x)

[Out]

-(84*a**6 + 560*a**5*b*x**2 + 1575*a**4*b**2*x**4 + 2400*a**3*b**3*x**6 + 2100*a**2*b**4*x**8 + 1008*a*b**5*x*
*10 + 210*b**6*x**12)/(1680*x**20)

________________________________________________________________________________________

Giac [A]  time = 1.14007, size = 95, normalized size = 1.13 \begin{align*} -\frac{210 \, b^{6} x^{12} + 1008 \, a b^{5} x^{10} + 2100 \, a^{2} b^{4} x^{8} + 2400 \, a^{3} b^{3} x^{6} + 1575 \, a^{4} b^{2} x^{4} + 560 \, a^{5} b x^{2} + 84 \, a^{6}}{1680 \, x^{20}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((b^2*x^4+2*a*b*x^2+a^2)^3/x^21,x, algorithm="giac")

[Out]

-1/1680*(210*b^6*x^12 + 1008*a*b^5*x^10 + 2100*a^2*b^4*x^8 + 2400*a^3*b^3*x^6 + 1575*a^4*b^2*x^4 + 560*a^5*b*x
^2 + 84*a^6)/x^20